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Recent experiment on a honeycomb-lattice Heisenberg antiferromagnet (AF) 
Bi3Mii40i2(N03) revealed a novel spin-liquid-like behavior down to low temperature, 
which was ascribed to the frustration effect due to the competition between the AF 
nearest- and next-nearest-neighbor interactions J\ and Ji. Motivated by the experiment, 
we study the ordering of the J\ -J% frustrated classical Heisenberg AF on a honeycomb 
lattice both by a low-temperature expansion and a Monte Carlo simulation. The model has 
been known to possess a massive degeneracy of the ground state, which, however, might 
be lifted due to thermal fluctuations leading to a unique ordered state, the effect known 
as 'order-by-disorder'. We find that the model exhibits an intriguing ordering behavior, 
particularly near the AF phase boundary. The energy scale of the order-by-disorder is 
suppressed there down to extremely low temperatures, giving rise to exotic spin-liquid 
states like a "ring-liquid" or a "pancake-liquid" state accompanied by the characteristic 
spin structure factor and the field-induced antiferromagnetism. We argue that the recent 
experimental data are explicable if the system is in such exotic spin-liquid states. 
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1. Introduction 

Recently, growing interest arises in the ordering of geometrically frustrated magnets. Par- 
ticular attention has been paid to possible 'spin-liquid' states stabilized due to the frustration 
effect, where spins remain to be disordered down to very low temperature without show- 
ing the standard magnetic long-range order. ^ Usually, geometrical frustration is realized in 
triangle-based lattices like the triangular, kagome and pyrochlore lattices combined with the 
antiferromagnetic (AF) coupling. By contrast, the honeycomb lattice, a hexagon-based lat- 
tice shown in Fig.l, is bipartite and is usually regarded as an unfrustrated lattice since it 
can accommodate the standard AF 'up-down' order. However, if the AF interaction works 
between the next-nearest-neighbor (nnn) sites in addition to the nearest-neighbor (nn) sites, 
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Fig. 1. (Color online) A honeycomb lattice, where J\ and Ji represent the nearest-neighbor (nn) and 
the next-nearest-neighbor (nnn) interactions, and a is the lattice spacing of the triangular Bravais 
lattice equal to the nnn distance of the honeycomb lattice. The green (red) arrows represent six 
equivalent nn (nnn) directions. 

the frustration effect arises due to the competition between the nn and the nnn couplings J\ 
and J2: See Fig.l. Since the honeycomb lattice is a loosely-coupled lattice with the number of 
the nn sites only three, it might be susceptible to the fluctuation effect caused by frustration, 
and its ordering property is of special interest. Quantum magnetism on the honeycomb lattice 
has attracted much attention. 2 ~ 6 ) Interest in the honeycomb lattice is also further promoted 
by the recent upsurging research interest in graphenes. 7 ' 8 ) 

An earlier theoretical study by Katsura et al on the honeycomb-lattice classical Heisenberg 
model with the nn and the nnn AF couplings ( J1-J2 model) revealed that, when the nnn AF 
coupling is moderately strong J2/J1 > 1/6, the classical ground state is infinitely degenerate, 
i.e., the ground-state manifold of the model consists of a set of spiral states characterized 
by generally incommensurate wavevectors q which form the ring surrounding the AF point 
in the wavevector space, rather than the discrete points. 9 ) If the nnn coupling is weaker 
J2/J1 < 1/6, on the other hand, the energy minimum occurs at the AF point. Such features 
are demonstrated in the Fourier-transformed energy of the model shown in Fig. 2(a) and (b). 
Based on this observation, Katsura et al suggested that the ground state of the model might 
be disordered at J2/J1 > 1/6 because the system might continue to fluctuate over these 
degenerate states characterized by different q. 

Interest in the honeycomb AF is promoted by the recent experiment by Azuma et al 
who observed that the S = 3/2 honeycomb-lattice Heisenberg AF Bi 3 Mn 4 0i2(N0 3 ) (BMNO) 
exhibited a spin-liquid-like behavior down to low temperature 0.4K, much lower than the 
Curie- Weiss temperature Tqw — ~ 257K. 10 ~ 12 ) The frustration effect due to the nnn coupling 
was invoked to explain the observed spin-liquid-like behavior. It thus remains most interesting 
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(a) 

J2/J1 = 0.18 




Fig. 2. (Color online) The minus of the Fourier-transformed energy of the model is plotted in the (h, I) 
plane, where q = —(h, I) is the wavevector, for the case of J2/J1 = 0.18 (a), and of J2/J1 = 0.30 
(b). The green ring represents the degenerate line of the ground state, which appears surrounding 
the AF point (0,2/ 

to understand the ordering process of the frustrated honeycomb-lattice Heisenberg AF. 

In fact, a naive expectation linking the observed infinite degeneracy to the disordered 
ground state needs careful examination, since such degeneracy is often lifted by fluctuations, 
either thermal or quantum, a phenomenon known as 'order- by-disorder', 13 ) which is often ob- 
served in frustrated spin systems. 14 ~ 16 ) In the honeycomb-lattice Heisenberg AF, one naturally 
expects that an infinite ring-like degeneracy of the ground state might be lifted by fluctua- 
tions, leading to a unique spiral ordered state characterized by a unique q. The candidate 
directions selected might be either the nn or the nnn directions as shown in Fig.l. 

We find that such an order-by-disorder mechanism is certainly operative in the honeycomb- 
lattice Heisenberg AF, leading to a unique ordered state where the threefold [Z3] directional 
symmetry of the honeycomb lattice is spontaneously broken. Usually, the order-by-disorder 
occurs at the energy scale of the main coupling, i.e., at the order of the Curie- Weiss tem- 
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perature Tew or the nn coupling J±, if somewhat suppressed due to the frustration effect. 
Therefore, the transition temperature to a unique ordered state selected by the order-by- 
disorder is not extremely low. An example might be a field-induced plateau phase realized in 
the Heisenberg AF on the triangular and kagome lattices. 17 ' 18 ) Emergence of the spin-liquid 
naively expected from the ground-state degeneracy might thereby be hampered in reality. 
Such limitations might particularly be severe in three-dimensional (3D) frustrated magnets. 

In such circumstances, an interesting possibility arises in the 2D honeycomb-lattice Heisen- 
berg AF. We observe that, in this model near its AF phase boundary J2/J1 ~ 1/6, the energy 
scale associated with the order-by-disorder is determined by J2 — | J\ rather than by J\ or J2 
itself, becoming arbitrary small around the AF phase boundary. Indeed, it tunes out that this 
regime, with the help of the ring-like degeneracy, gives rise to a variety of interesting behaviors 
such as exotic types of spin-liquid states, which we call "ring-liquid" and the "pancake-liquid" 
states, as well as the field-induced AF order. 

2. Model and method 

The model we consider is the classical Heisenberg AF on the 2D honeycomb lattice, whose 
Hamiltonian is given by 

n = Ji^2s i -s j + j 2 ^2s i -s j , (i) 

{{id)) 

where the sum is taken over all nn and nnn pairs for J\ and J2, respectively. We study the 
ordering properties of the model both by a low-temperature expansion and by a Monte Carlo 
(MC) simulation. 

The low-temperature expansion is made following the method of Ref., 16 ) expanding around 
an arbitrary state in the ground-state manifold at the harmonic order. Some of the details is 
given in Appendix. 

MC is performed on the basis of the standard heat-bath method combined with the over- 
relaxation technique and the temperature-exchange method. Various values of j'2 = J2/J1 in 
the range j'2 = [0.17, 0.5] are studied, employing both periodic and free boundary conditions 
(BCs). The lattice contains 2 x L 2 spins, with L ranging between 24 < L < 72. In the present 
paper, we study mainly the case of J2 = 0.1830 • • • , which corresponds to l/[4cos | + 2] with 
1 = T§ ' f° r s y s t ems under periodic BCs. This particular value of j'2 is chosen to minimize 
the finite-size effect due to the mismatch between the incommensurability of the helix and 
the applied periodic BCs: The helix along the nnn direction becomes commensurate with the 
lattice at this value of j'2 if the L- value is chosen as multiples of twelve. 
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Fig. 3. (Color online) The J2/ Ji-dependence of the free-energy difference between the spiral states 
running along the nn direction and along the nnn direction, AF = i^nnn) — F(nn), calculated by 
the low-temperature expansion. At Jij J\ 5 0.232 (c 0.232), the nnn (nn) direction has the lower 
free energy, whereas the difference is small for J2/J1 ~ 0.232. 



3. Results 

First, we report on our results of the low-temperature expansion. We find that thermal 
fluctuations select a particular incommensurate spiral state whose helix axis runs along the 
nnn direction for J2/J1 = ji 5 0.232, but along the nn direction for j'2 <i 0.232: See Fig. 3. 
The calculated free-energy difference between these two directions, AF, is shown in Fig. 3 as 
a function of J2/J1. Over the range 1/6 < j'2 5 0.232, AF stays small, suggesting that the 
selection is weak, whereas, for j'2 0.232, AF is large. Our result of thermal selection differs 
from that of the recent T = calculation of Ref. 19 ^ for the quantum honeycomb J\-J% model, 
where the nn direction was selected irrespective of the jVvalue for 1/6 < j'2 < 0.5. 

Next, we present our MC results. In Fig. 4(a), we show the specific heat per spin for 
22 = 0.1830. (Here and below, the unit of the energy and the temperature is taken to be 
J\.) The specific heat exhibits three peaks, a broad peak around T ~ 0.12 associated with 
the onset of the short-range order, a sharp diverging peak at T = T c ~ 0.013 associated 
with a phase transition, and another peak at a lower temperature T = T* whose position 
and height depend on the system size L considerably. The transition at T c seems continuous. 
The Curie- Weiss temperature is estimated to be Tew — —1-68: See the inset of Fig.4(b). The 
diverging peak at T c ~ 0.013 is associated with the spontaneous breaking of the Z3 directional 
symmetry of the lattice. 

This can be confirmed from the order parameter 7713 describing the Z3 directional- 
symmetry breaking. This order parameter is defined by 

m 3 = (|m 3 |), ^3 = eiei + e 2 e 2 + e 3 e 3 , (2) 
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Fig. 4. (Color online) The temperature and size dependence of several physical quantities at J%j J\ = 
0.1830. (a) The specific heat per spin. The inset represents a wider temperature range, (b) The 
threefold directional-symmetry-breaking order parameter m^. The inset represents the Curie- Weiss 
plot, i.e., the temperature dependence of the inverse susceptibility per spin. The Curie- Weiss 
temperature is estimated to be Tew — —1-68. 



where e\ = (0,1), £2 = ( — 2 ' — 2) an< ^ ^ 3 = (2 > — 2)' anc ^ e M 0" = 1)2,3) is the total nn 
bond energy (normalized per bond) along the direction. When the spiral runs along the 
nnn direction, 7713 takes a value 1 — cos | at T = 0. For J2/J1 = 0.1830, the saturation value 
is then m^(T = 0) = 1 — ~ 0.134. As can be seen from Fig. 4(b), the directional order 
sets in at T c ~ 0.013 which coincides with the specific-heat peak. Even below T c , 1113 does not 
saturate but increases further, and saturates below T* . 

In Figs. 5 and 6, we show the spin structure factor F(q) in the (h, I) plane for the case 
of J2/J1 = 0.1830, where q = %(h,l) and a is the lattice constant of the underlying Bra- 
vais lattice: See Fig.l. The spin structure factor is defined by F(q) = (\S q \ 2 ^, where (• • • ) 
represents the thermal average and S q is the Fourier transform of the spins. Fig. 5 represents 
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F(q) at a temperature slightly above T c shown in the entire q space around the V point, while 
Fig. 6 is a magnified view of a part of the q space around an AF point (0,2/\/3) at several 
temperatures across T c and T* . At a temperature below T*, F(q) exhibits sharp point-like 
spots in the nnn direction, consistently with the prediction of the low-temperature expansion. 
In fact, the nnn direction is chosen at lower temperatures for 1/6 IS 32 ~ 0.232 for large 
enough lattices under periodic boundary conditions. However, the other nn direction tends to 
be stabilized for smaller sizes. Finite-size effects seem significant here, presumably reflecting 
the small free-energy difference between the two directions shown in Fig. 3. 



T=0.0150 




-1.5 -1 -0.5 0.5 1 1.5 



h 

Fig. 5. (Color online) The intensity plot of the spin structure factor in the entire (h, I) plane around 
the r point (0,0) for the case of J2/J1 — 0.1830, where q — —(h,l). The temperatures is T = 
0.0150 and the lattice size is L — 72. In this simulation, we turn off the temperature-exchange 
process when measuring F(q) to probe the way of the symmetry breaking clearly, while it is turned 
on in equilibrating the system prior to the measurement. 

At a temperature between T* and T c , F(q) still selects the state in the nnn direction, but 
now the intensity becomes diffuse exhibiting a ridge along the ring direction: See Fig. 6(b). It 
means that the fluctuation in the q- vector direction becomes enhanced. At temperatures above 
T c , the Z3 symmetry is restored and F{q) exhibits a pronounced ring-like pattern surrounding 
the AF point, which faithfully reflects the ground-state degeneracy: See Fig. 6(c). On further 
increasing the temperature, the center of the ring is 'buried', giving rise to a "pancake-like" 
pattern shown in Fig. 6(d). Note that the radius of such a "pancake" is still determined by the 
radius of the degenerate ring, and hence, hardly depends on the temperature. 

In Fig. 7(a), we show for 32 = 0.1830 a segment of F(q) along the h = axis at various 
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Fig. 6. (Color online) The intensity plot of the spin structure factor in the (h,l) plane for J2/J1 = 
0.1830, magnified for a part of the plane around an AF point (0,2/^/5). The temperatures are 
T = 0.0070 (a), T = 0.0090 (b), T = 0.0150 (c), and T = 0.0694 (d), respectively. The lattice 
size is L — 72. We turn off the temperature-exchange process when measuring F(q) to probe 
the way of the symmetry breaking clearly, while it is turned on in equilibrating the system prior 
to the measurement. The next-nearest-neighbor (nearest- neighbor) directions of the honeycomb 
lattice are given by the green solid (broken) lines in (a) and (b). The ring corresponding to the 
ground-state degeneracy is shown by the broken circles in (c) and (d). 



temperatures. F{q) exhibits the double-peak pattern corresponding to the ring for 0.013 
T "S 0.06 and the flat pancake-like pattern for 0.06 "S T 'S 0.11. Note that the pancake 
pattern is not fittable by the double-Lorenzian centered at the ring positions. These states 
with characteristic shapes of F{q) are all paramagnetic states in the sense that no symmetry 
breaking occurs there, but they differ from the standard paramagnetic state in the sense that 



8/17 



J. Phys. Soc. Jpn. 



Full Paper 




Fig. 7. (Color online) The spin structure factor in the "ring-liquid" and the "pancake-liquid" regions, 
(a) The i-dependence of the spin structure factor along the h = line for various temperatures 
above T c at J2/J1 = 0.1830. The lattice size is L = 72. The "ring-liquid" state is realized at 
0.013 < T < 0.06, while the "pancake-liquid" state is realized at 0.06 < T < 0.11. The inset 
represents the temperature dependence of the width (at half maximum) of the -F^cd^D-peak 
for several J2/ Ji-values. The crossover points between the ring-liquid and the pancake-liquid, 
and the ones between the pancake-liquid and the standard paramagnet are given by arrows, (b) 
The |q|-dependence of the powder-averaged spin structure factor, i.e., the spin structure factor 
averaged over angles with the modulus |q|, in the ring- liquid and the pancake- liquid regions for 
J2/J1 = 0.17 and 0.1830. The lattice size is L — 72. The inset exhibits the J2/«/i-dependence of 
the 'antiferromagnetic correlation length' £af m units of a corresponding to the inverse width of 
the F ave (|q|)-peak, which is given by the radius of the degenerate ring. 



the associated F(q) exhibits quite unusual shape. We call these exotic paramagnetic states a 
"ring-liquid" and a "pancake-liquid" state. 

In principle, it should be possible to experimentally observe such characteristic changes 
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of F(q) by means of neutron scattering. Since only powder samples have been available so 
far for BMNO, we also calculate the powder averaged F{q) for several typical cases, and 
the results are shown in Fig.7(b). The power-averaged F(q) exhibits a peak centered around 
the AF |q|-point with the width equal to the radius of the degenerate ring of the ground- 
state manifold. Our present F(q) seems fully consistent with the recent experimental data 
on BMNO. 11 ) In our model calculation, the width of the peak in the "ring-liquid" or the 
"pancake-liquid" regimes, which would be interpreted as the inverse 'AF correlation length' 
£af experimentally, remains essentially temperature- independent. This is demonstrated in 
the inset of Fig.7(a) where the temperature dependence of the width of the F ave (|g|)-peak 
is shown for several J2/ J\- values. In the ring and the pancake regimes, the width exhibits 
negligible or even opposite temperature dependence. In the inset of Fig. 7(b), we show the 
■h/ J\ -dependence of the 'AF correlation length' ^af '■ Towards the AF boundary 22c = 1/6, 
£af diverges as oc \j\J 22 — 22c- Nevertheless, it stays short even fairly close to the AF phase 
boundary: For example, £af — a at 22 = 0.1830, and £af — 2a at 22 = 0.17. The recent 
experiment on BMNO yields £af — 8A ~ 1.6a. 11 ) 

Now, we touch upon the issue of T*, a transition-like temperature manifesting itself in 
the lower specific-heat peak and in the sharp rise of the directional order parameter m%. We 
deduce that T* corresponds to the Z2-vortex binding-unbinding transition. It has been known 
that the frustration-induced noncollinear order of the Heisenberg spin sustains a characteris- 
tic topological excitation called a Z2 vortex. 20 ) Indeed, recent studies have indicated that a 
topological transition driven by the binding-unbinding of such Z2 vortices occurs in certain 
frustrated 2D Heisenberg models with the noncollinear spin order. 21 ' 22 ) Our preliminary MC 
study of the Z2-vortex distribution supports such an identification. Further details will be 
reported elsewhere. 

Next, we study the effect of applied magnetic fields. In the vicinity of the AF phase 
boundary, fields tend to induce the AF order. This is demonstrated in Fig. 8(a) where the 
AF order parameter vtiaf is plotted versus the field intensity for several temperatures at 
J2/J1 = 0.175 close to the AF phase boundary. The AF order parameter is defined by 

m A F = l(\(m x A y-m^)\), (3) 

where m A y and m X g are the transverse components (perpendicular to the applied magnetic 
field) of the sublattice magnetization per spin of the sublattices A and B. One can see that 
the AF correlation is significantly enhanced by fields in the temperature range above T c . Note 
that rriAF should vanish in the thermodynamic limit in the 2D Heisenberg model, yet gives a 
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Fig. 8. (Color online) (a) The magnetic field dependence of the AF order parameter tyiaf- (b) The 
field dependence of the spin structure factor of the transverse components versus I with h — for 
J2/J1 = 0.175. The temperature is T = 0.042 above T c . The lattice size is L = 36. Applied fields 
induce a sharp peak around the AF point, I = 2/\/3. 

measure of the AF correlations in finite systems. In real systems, a weak 3D coupling might 
eventually realize the true AF long-range order. The field dependence of the spin structure 
factor is also shown in Fig. 8(b). The region of such a field-induced AF is also depicted in 
Fig. 9(b), which occupies the region around the pancake-liquid regime. It would be natural to 
expect that the pancake- liquid, which may be regarded as a modified ring-liquid state where 
the AF component is enforced, is favored in the AF-state formation. We note that a similar 
field-induced AF has recently been reported on BMNO by Matsuda et al. 11 ^ 

Our results are summarized in the J2/J1 versus temperature phase diagram of Fig. 9(a). 
In Fig. 9(b), the vicinity of the AF boundary is magnified by extending the temperature 
range to higher temperatures. Toward the AF phase boundary, the transition line tends to 
T = 0, stabilizing the spin-liquid state without any symmetry breaking down to very low 
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temperature. Indeed, T c /\Tcw\ is estimated to be as low as ~ 0.008 at 32 = 0.1830 and 
~ 0.001 at j2 = 0.17. In the vicinity of the AF phase boundary, the paramagnetic state is 
of unusual type characterized by the ring-like or the pancake-like spin structure factors, the 
"ring-liquid" or the "pancake-liquid" state. 

It should be emphasized that the strong frustration effect and the subsequent massive 
degeneracy alone are not enough to stabilize such spin-liquid states down to low temperature, 
since the order-by-disorder mechanism usually comes into play leading to the selection of 
a particular ordered state at the energy scale of the main exchange interaction 0(J\). In 
the vicinity of the AF phase boundary of this model, the AF state, which remains to be a 
locally unstable state in itself, conspires with the frustration-induced massive degeneracy to 
stabilize the special types of spin-liquid states down to extremely low temperature. The ring- 
like degeneracy in the q-space also plays an essential role in realizing novel spin-liquid states 
here, which is absent in, e.g., the J1-J2 Heisenberg model on a square lattice, another 2D spin 
system exhibiting the order-by-disorder phenomena. 23 ) 

If the J2/Ji-value is increased away from the AF phase boundary, T c increases, whereas 
the order of the transition changes from continuous to first-order at around 32 ~ 0.25. It 
remains first-order for 0.25 5 32 ~ 0.35 and becomes continuous again beyond 32 ~ 0.35: See 
Fig.9(a). Similar first-order Z3 directional-symmetry-breaking transition was reported also for 
certain triangular-lattice Heisenberg AF. 24 ) The behavior of the T c -line is consistent with the 
recent calculation by Mulder et al, 19 ^ whereas the change in the order of the transition or the 
existence of another T*-line was not reported there. 

4. Discussion and summary 

We wish to emphasize that our results are consistent with the recent experiment on 
BMNO 10 ~ 12) if BMNO lies close to the AF phase boundary, (i) The spin-liquid-like state 
is stabilized down to very low temperature. Our model provides a mechanism to significantly 
suppress the ordering in the vicinity of the AF phase boundary, (ii) Neutron-scattering inten- 
sity observed for powder samples, which exhibits an extremely broad peak centered around 
the AF point with apparently short AF correlation length, is consistent with the correspond- 
ing pattern expected for the pancake-liquid (or the ring-liquid) state of the present model, 
(iii) Experimentally observed field-induced AF is fully consistent with the corresponding phe- 
nomena observed in the pancake-liquid (or in the ring-liquid) state of the present model, (iv) 
Overall behavior of the specific heat and the susceptibility above T c are consistent with the 
experimental results. 10 ) 
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9. (Color online) Phase diagram of the model, (a) The T c -curvc representing the Z3 directional- 
symmetry-breaking transition is plotted versus J2/J1, together with the T*-curve determined 
from the lowest specific-heat-peak position for the largest size studied. Note that T* still depends 
on the system size appreciably so that the precise location of the T*-curve is uncertain. We 
deduce that T* corresponds to the ^-vortex binding- unbinding transition (see the text). The 
directional-symmetry-breaking transition at T c is continuous in the smaller and in the larger 
J2/J1 regimes, but is first-order in the intermediate J2/J1 regime. There is a tendency that the 
directional-symmetry-breaking transition, which is originally continuous, becomes of first-order 
when it merges with the T*-line. (b) represents the vicinity of the AF phase boundary J2/J1 = g 
including the higher temperature range above T c . The regions of the "pancake-liquid" , the "ring- 
liquid" and the standard paramagnetic states are indicated, each of which is separated by a 
crossover line. The crossover line between the ring-liquid and the pancake liquid is determined 
from the criterion whether F(q) exhibits a single peak or double peaks, while that between the 
ring-liquid and the standard paramagnet is determined from the temperature dependence of the 
width of the i 7, avo (q)-peak given in the inset of Fig. 7(a). The region exhibiting the field-induced 
AF is indicated by shaded area, which is determined from the criterion that the maximum value 
of rriAF induced by applied fields is greater than that in zero field by more than 50%. 
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Of course, in real experiments, several points not taken into account in the present model 
might play a role. Here we discuss the following two points, (a) Real BMNO is a bilayer 
honeycomb system with the AF interlayer coupling. We have also performed a preliminary 
simulation on the bilayer model, to find that the main results, particularly the points (i) ~ 
(iv) above, do not change, at least qualitatively, (b) Real samples are likely to possess certain 
amounts of defects and impurities which might cause the glassy behavior at low temperatures, 
masking an intrinsic behavior of the pure system. For example, the expected phase transition 
has not been observed so far for BMNO down to 0.4K. The reason of this might be either (i) 
actual T c is lower than 0.4K, or (ii) the system at low temperatures is gradually stuck into 
the glassy state because of the randomness. In any case, further measurements on purer, and 
possibly, single crystal is desirable to clarify the issue. 

In summary, we have revealed an intriguing ordering behavior of the frustrated honeycomb- 
lattice Heisenberg AF. Near the AF phase boundary, exotic spin-liquid states like the ring- 
liquid or the pancake- liquid state are stabilized down to extremely low temperatures, accom- 
panied by the characteristic spin structure factor and the field-induced antiferromagnetism. 
Our results seem consistent with the recent experimental data on BMNO. 
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Appendix: Low-temperature expansion 

Here we calculate the partition function Z of the model described by the Hamiltonian %, 



by means of a low-temperature expansion from an arbitrary state in the ground-state manifold. 
Let the ground-state spin orientation be Si, which is assumed to lie in the xy plane. Si is an 
incommensurate spiral state characterized by the wavevector q belonging to the degenerate 
ring of the ground-state manifold. By introducing the deviation vector 7Tj, which satisfies 
7Tj _L Si, one has 



Z 




(A-l) 




By decomposing the 7Tj-vector into the xy- and z-components, 
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■Ki = e z <t>i + [e z x Si]xi 



and expanding the Hamiltonian up to the quadratic order both in 7T and (f>, the partition 
function Z 



Z = J D(pD X e~ m 



(A-2) 



can be evaluated by the Gaussian integrals. Neglecting the terms independent of the wavevec- 
tor q specifying the ground state, which even include the divergent terms, we finally get the 
following expression of the g-dependent part of the free energy, 



F(q)/T = J dg / (In[Wii(g,g / ) + |Wi2(9,9 / )|] 



+ ln[W 11 (q,q')-\W 12 (q,q')\}). (A-3) 



where 



W 11 (q,q')=2J 2 



cos q x cos q x 



Qx + V3q' y q x + \/3q y 
+ cos cos - 



Qx ~ V3q' y q x - y/3q y 
+ cos cos 



X(q), (A-4) 



W l2 (q,q') = J l 



cosa(g)e 



.(l' X +1y/V3) 



+ cos (a(q) - ( k±^ )e </^ 



+ cos(a(q) - q x )e % 2 



The function \(q) is defined by 



with 



\{q) = 2J 2 e{q) - J 1 ^3 + 2e{q), 



, , , Qx + V3q y , Qx ~ V3q y 
e(q) = cos q x + cos + cos ■ 



(A-5) 



(A-6) 



(A-7) 



2 2 

while a(g), representing a phase angle between the two nn spins in a unit cell belonging to 
the two sublattices A and B as o\ A ^ = q ■ r\ A ^ and = q • r\ A ^ + a(q), is given by 



cos a(q) = — (1 + cos q x + cos 



q -^^)/^THq), 



(A-8) 
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sma(q) = - (sin ga + sin v Hy )/y/3 + 2e(q). (A-9) 
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